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Abstract. We investigate the ionization structure in the non-spherical winds 
and disks of B[e] stars. Especially the luminous B[e] supergiants seem to have 
outflowing disks which are neutral in hydrogen already close to the stellar surface. 
The existence of neutral material so close to the central star is surprising and 
needs to be investigated in detail. We perform our model calculations mainly in 
the equatorial plane, trying to find a plausible scenario that leads to recombina- 
tion in the vicinity of the hot stars. Two different approaches are presented that 
both result in a hydrogen neutral equatorial region. We especially focus on the 
influence of stellar rotation which is known to play a significant role in shaping 
non-spherical winds. We show that a rotating star can have a neutral equatorial 
wind, even without the need of a density enhancement due to bi-stability or wind 
compression, but simply due to graviational darkening in combination with the 
known high mass- loss rates of the B[e] stars and especially the B[e] supergiants. 



1. Introduction 



The circumstellar material of B[e] stars is known to be non-spherically symmet- 
ric, and the best studied sample of B[e] stars, the B[e] supergiants, even show 
strong indication for the existence of an equatorial gas and dust disk. The hy- 
brid character of their optical spectra have been interpreted by a non-spherical 
wind with a typical line-driven wind in p olar, and a slow, high-density wind in 
equatorial direction ( Zickgraf et alJ[T"985l . and this volume). Indications for the 
existence of a neutral gas a nd dust disk come (i) fr om molecular emission like 
the CO first overtone bands ( McGregor et al.lll988al|b|h (n) the e xistence of dust 
causing a strong mid-IR excess emission (see e.g. Ivborsl I1999L and references 
therein), and (iii) the detection of strong [Oi] emi ssion in the high-resolution 
optical spectra of severa l B[e] stars and supergiants (|Kraus fc Borges Fernandesl 
2005fc iKraus et al 1 l2005l . and poster presentation, this volume), the latter even 
speaks in favour of a huge amount of hydrogen neutral material in the vicin- 
ity of these luminous stars. This neutral material can only exist in the dense 
disks around these luminous stars, and model calculations have shown that these 
disks must be neutral already near the stellar surface to account for the observed 
line luminosities and to keep the needed disk mass loss rate at a reliable value 
((Kraus Borges Fernandesl l2005 i, and Kraus et al. poster presentation, this 
volume) . The goal of our study is therefore to find scenarios which allow neutral 
material to exist close to the surfaces of these stars. 
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Figure 1. Distribution of the mass flux (left), escape (and hence terminal, 
because Voo ~ v CS c) velocity (mid) and hydrogen density (right) over the 
stellar surface. The parameters are normalized to their polar values. Different 
line st yles are for different model parameters as specified in lKraus fc LamersI 
(200]|. 




2. Ionization structure calculations in non-spherical winds 

We present two different approaches which both lead to the formation of a 
neutral equatorial region. The first one is an empirical model which assumes a 
non-spherical wind with a surface mass flux that increases from pole to equator. 
The second one deals with the influence of rigid stellar rotation. 

2.1. Empirical non-spherical wind model 

In a first attempt we dev elopped a code to calcu late the ionization structure in 
a pure H plus He wind ( Kraus k, LamersI 120031 ) . We used empirically defined 



latitude dependent mass flux and escape (or terminal) velocity distributions 
(left and mid panel of Fig. 1), which define a surface density distribution that 
accounts for the observed velocity and density ratios between equator and poles, 
and creates a disk-like structure with an opening angle of 20° — 30° (mid and 
right panel of Fig. 1). 

We solved the ionization balance equations for H and He in the on-the-spot 
approximation which means that all photons generated via recombination and 
suitable for re-ionization are absorbed in situ. The electron temperature and ra- 
dial wind velocity were set constant and we used the point source approximation 
for the stellar radiation field. These restrictions resulted in the calculation of an 
upper limit of the recombination distances which are plotted for three different 
model parameters in Fig. 2 (for more model details see Kraus & Lamers 2003). 

All these models for which the total mass loss rates were in the range 2.4 x 
10~ 6 ...2.4 x lO" 5 M yr _1 , which is a reasonable range for B[e] supergiants, 
resulted in a hydrogen neutral disk close to the hot stellar surface. 

2.2. Stellar rotation and the formation of a non-spherical wind 

During the last years, the non-negligible influe nce of stellar rotation on stellar 
evolution has been investigated in detail (see e.g. lMaederll999l : lMaeder k, Mevnetl 
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Figure 2. Three examples for the ionization structure of H (dashe d) and 
He (solid) in non-spherical winds (taken from lKraus fc Lamersll20 03). The 
material to the right of the full or dashed line is neutral. Both elements 
clearly recombine close to the stellar surface in the equatorial plane, leading 
to a hydrogen neutral disk. 

2000 ; iMevnet fc Maederll2000l . and this volume) , and it has been suggested that 
the appeara nce of non-spherical winds around lumino us stars might be caused 
by rotation (|Maederl 120021 : iMaeder k, Desiaauesl l2001). We therefore started to 
investigate the influence of rigid stellar rotation on the formation of a hydrogen 
neutral disk. Rotation causes a flattening of the stellar surface and therefore 
a reduction of the local net gravity in the equatorial region. The decrease in 
gravity from pole to equator equally results in a decrease of the stellar flux 
which is proportional to the local net gravity. Hence, the effective temperature 
also decreases from pole to equ ator (left panel o f Fig.|3J), known as the gravity 
darkening (or polar brightening, von Zeine]lll924l ). 

The latitude dependence of the gravity and effective temperature has also 
impacts on the stellar wind parameters: The escape velocity following from the 
balance between gravitational and centripetal forces becomes latitude depen- 
dent, decreasing from pole to equator. The same holds for the terminal wind 
velocity which is (for line-driven winds) proportional to the escape velocity. And 
even the mass flux from the star tends to decrease from pole to equator when 
implem enting; gravity darkening into the CAK theory, as shown bv lOwocki et al.l 
( 1998ft . 
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Figure 3. Normalized surface distributions of the effective temperature 
(left) and the wind density (right) as a function of the co-latitute, 0. The 
curves from top to bottom are for an increasing rotational velocity indicated 
by the parameter ui, which is the ratio of the equatorial rotation velocity over 
the critical velocity. 



More important for the ionization structure calculations is the density in the 
wind. According to the mass continuity equation the hydrogen density scales as 
nu ~ F m /(vr 2 ). A proper elaboration of this ratio, taking into account the ro- 
tational distortion of the stellar surface, resulted in surface density distributions 
of rotating stars as shown in Fig. |3 (right panel). From this plot it is clear, that 
the density at any given distance also decreases from pole to equator. A rotating 
star, when neglecting additional effects like bi-stability or wind compression, will 
therefore have a less dense wind in the equatorial region. 

Both important parameters in the ionization balance equations, i.e. the ef- 
fective surface temperature and the surface density, decrease from pole to equa- 
tor. While the decrease in surface temperature tends to decrease the number of 
available photons suitable to ionize H and He, the decrease in surface density 
reduces the optical depth along the line of sight from a point in the wind to the 
star. Both effects are therefore counteracting with respect to the location where 
recombination takes place: While a reduction of ionizing photons will shift the 
recombination distance towards the star, the reduction in optical depth along 
the direction to the star will shift it further outwards. The outcome of the 
ionization balance equations is therefore unpredictible and very sensitive to the 
choosen input parameters. This makes it worth to investigate the ionization 
structure in the wind of a rotating star in more detail. 

In our model calculations, the wind is assumed to be composed of H and 
He, and we solve the ionization balance equations in the equatorial plane, only. 
The radiation field from the star is calculated properly in 2D, which means that 
for the available radiation at every point in the equatorial wind we integrate 
the stellar radiation over the non-spherical surface and take into account the 
latitude variations of the individual parameters like effective temperature, net 
gravity, escape velocity and mass flux. We make, however, several assumptions 
and simplifications for the calculations which all have in common that they 
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Figure 4. Left: Stellar luminosity as a function of u) (w = u r ot/^crit) for 
our model stars with T e ff iP oic = 17 000 K and i?„ = 82i? Q . Right: The ratio 
between the mass loss rate, integrated over the stellar surface, and the polar 
mass flux as a function of u>. The decrease in both, luminosity and total mass 
loss rate, is about a factor of 2 between the non-rotating and the critically 
rotating star. 



shift the location of recombination to larger distances. This means that we are 
calculating an upper limit for the recombination distances 1 . 

We performed a series of test calculations where we fixed the polar values 
of the stars, i.e. the polar effective temperature and gravity which determine 
the radiation temperature of the ionizing stellar radiation, and the polar ter- 
minal velocity. Fixing the polar effective temperature means that we are not 
calculating a star that is spinning up which would result in an increase of polar 
temperature with increasing rotation velocity, but we are calculating stars with 
the same polar effective temperature having different rotation velocities. This 
means that we are dealing with stars of different luminosities. 

The difference in luminosity of these stars can be seen in Fig.|I] (left panel) 
where we plotted the distribution of the stellar luminosities of our model stars as 
a function of stellar rotation. This difference is largest between the non-rotating 
and the critically rotating star, and is about a factor 2. 

Also shown in Fig.Q] (right panel) is the variation of the total mass loss 
of the rotating stars. For the choosen polar mass fluxes of 1.5 x 10 -6 and 
1.5 x 10 _5 gs _1 cm -2 , the mass loss rates decrease from 10~ 5 and 10~ 4 MQyr _1 
(non-rotating) to 5 x 10 -6 and 5 x lO _5 M0yr _1 (critically rotating), respectively. 

First results of the ionization structure calculations in the equatorial winds 
of luminous rotating stars are shown in Fig.|SJ These calculations have been 
performed for polar mass fluxes that differ by a factor 10 to show clearly the 
influence of the input parameters. For the high mass flux (bottom panels) both 
elements recombine right above the stellar surface for to > 0.8, resulting in a 
(hydrogen) neutral equatorial outflowing wind. The lower mass flux model (top 
panels) also shows the trend of a decrease in recombination distance with in- 
creasing rotation. In the case of hydrogen, this distance has however a minimum 



1 For details on the assumptions and their influence see lKraus fc Lamersl (2005). 
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Figure 5. Helium (left) and hydrogen (right) recombination distances in 
the equatorial plane as a function of to. The polar mass flux used for each 
model is indicated. While for the high mass flux model hydrogen recombincs 
right above the stellar surface for u> > 0.8, the recombination distance for the 
model with lower mass flux has a minimum at about u) ~ 0.7 and increases 
again for higher rotation rates. 



around to ~ 0.7, and increases again for higher rotation rates. This behaviour is 
unexpected but can be understood by looking at the variation of the individual 
equatorial surface parameters, especially of the temperature and wind density, 
with increasing rotation velocity as shown in Fig-IHl This figure shows that with 
increasing rotation velocity, the density drops much quicker from pole to equator 
than the temperature. For recombination to take place right above the stellar 
surface, the number of ionizing photons has to be reduced. This can be done 
either by decreasing the radiation temperature, or by increasing the equatorial 
surface density and hence the optical depth. Since the decrease in radiation 
temperature is determined by the rotation velocity (with a fixed input polar 
value) we can only increase the input polar mass flux to achieve a higher surface 
density for a given rotation velocity. The density in the model shown in the top 
panel of Fig. [5] is not high enough any more for stars with to > 0.7 to absorb 
the ionizing photons provided by the still rather high radiation temperature. 
Therefore, ionization takes over again and shifts the recombination distance for 
higher rotation velocities further out. Whether recombination takes place close 
to the star therefore sensitively depends on the choosen input parameters of the 
rotating star. 
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Figure 6. Surface parameters in the equatorial plane as functions of the 
rotation velocity. «h is the density at any distance. 



3. Conclusions 

We presented two different wind scenarios that resulted in the formation of a 
hydrogen neutral region in the equatorial plane: a non-spherical wind model with 
a surface mass flux that increases from pole to equator, and a model investigating 
the influence of rigid rotation on the wind structure. We show that the model of 
a rotating luminous star can have a hydrogen neutral equatorial region simply 
due to the effects of gravitational darkening. In that case, it makes however no 
sense to speak of a "disk" because the density in the equatorial plane (when 
neglecting additional effects like bi-stability and/or wind compression) is much 
lower than in the polar wind. 

Acknowledgments. M.K. acknowledges financial support from the Ned- 
erlandse Organisatie voor Wetenschappelijk Onderzoek grant No. 614.000.310. 

References 

Kraus, M. & Borges Fernandes, M. 2005, in The Nature and Evolution of Disks Around 

Hot Stars, ed. R. Ignace & K.G. Gayley, (San Francisco: ASP), p. 254 
Kraus, M. & Lamers, H.J.G.L.M. 2003, A&A, 405, 165 
Kraus, M. & Lamers, H.J.G.L.M. 2005, submitted to A&A 

Kraus, M., Borges Fernandes, M., de Araiijo, F.X. & Lamers, H.J.G.L.M. 2005, A&A, 

441, 289 
Maeder, A. 1999, A&A 347, 185 
Maeder, A. 2002, A&A 392, 575 
Maeder, A. & Desjaques, V. 2001, A&A 372, L9 
Maeder, A. & Meynet, G. 2000, A&A 361, 159 
McGregor, P.J., Hyland, A.R. & Hillier, D.J. 1988a, ApJ 324, 1071 
McGregor, P.J., Hillier, D.J. & Hyland, A.R. 1988b, ApJ 334, 639 
Meynet, G. & Maeder, A. 2000, A&A 361, 101 
Owocki, S.P., Cranmer, S.R. & Gayley, K.G. 1998, Ap&SS 260, 149 
Voors, R.H.M. 1999, PhD thesis, Utrecht University 
von Zeipel, H. 1924, MNRAS 84, 665 

Zickgraf, F.-J., Wolf, B., Stahl, O., Leitherer, C. & Klare, G. 1985, A&A 143, 421 



